
STABILITY AND DRAG CHARACTERISTICS 
AT MACH NUMBERS OF 10 AND 26 OF A 
PROPOSED SLENDER ATMOSPHERIC PROBE 


by Gerald N. Malcolm 

Ames Research Center 
Moffett Field, Calif. 




NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


APRIL 1967 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



NASA TN D-3917 


STABILITY AND DRAG CHARACTERISTICS AT MACH NUMBERS OF 

10 AND 26 OF A PROPOSED SLENDER ATMOSPHERIC PROBE 

By Gerald N. Malcolm 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



STABILITY AND DRAG CHARACTERISTICS AT MACH NUMBERS OF 


10 AND 26 OF A PROPOSED SLENDER ATMOSPHERIC PROBE 
By Gerald N. Malcolm 
Ames Research Center 


SUMMARY 


Free -flight tests at Mach numbers 10 and 26 were conducted on models of 
a slender, flare -stabilized cone cylinder proposed by Goddard Space Flight 
Center for use as a Mars -atmosphere probe. The tests were conducted to deter- 
mine drag, stability, and flow-field characteristics in both nose-forward and 
nose-rearward attitudes. Results of nose-forward tests at Mach numbers of 10 
and 26 and at Reynolds numbers of 0.6xl0 6 showed extremely nonlinear but sta- 
ble pitching -moment characteristics. The stability at low angles of attack 
(due, at least in part, to the observed laminar -boundary-layer separation over 
the cylinder and flare) decreased rapidly with increasing angle of attack 
until the bow-shock wave approached the windward side of the flare. Further 
increasing the angle of attack produced a notable increase in stability. No 
significant Mach number effects were observed in either the stability or drag 
results. Tests results at a Mach number of 10 showed the configuration was 
unstable in the nose-rearward attitude. The pitching moment for the model in 
this attitude was deduced from the rotational acceleration of the model. The 
pitching moment for the entire angle -of -attack range from a = 0° to a = l80° 
was determined by extrapolation of the experimental results from nose -forward 
and nose -rearward tests plus an iterative improvement of the extrapolation by 
a technique which matches a synthesized motion to several cycles of an 
observed free -flight motion. Models of the probe were found to be dynamically 
stable at all angles of attack tested. 


INTRODUCTION 


The design of an instrumented soft lander for exploration of the Martian 
surface is strongly dependent upon the characteristics of the Martian atmos- 
phere through which it must fly and which decelerates the lander helping to 
bring it to rest at the surface. Atmospheric probe missions to precede the 
lander mission in which the entry vehicle would, in effect, take soundings of 
the Martian atmosphere and relay the measurements to Earth prior to impact are 
therefore being considered. 

One approach, proposed by Goddard Space Flight Center, would use a 
slender flare -stabilized entry vehicle on which measurements of impact and 
static pressures, accelerations, and mass spectra of gas samples would be 
interpreted in terms of atmosphere characteristics. In support of the Goddard 
effort, the hypersonic aerodynamics of models of the slender -body probe have 



been investigated at the Ames Research Center in the prototype hypervelocity 
free-flight facility and the pressurized ballistic range. Two points were of 
particular interest. The first was to verify that the probe configuration was 
statically unstable in the nose-rearward attitude and had only one stable 
nose -forward trim attitude. This is important since the body would be enter- 
ing the atmosphere passively and therefore could be at any angle of attack. 

The test conditions corresponded to a combination of the highest Mach number 
and lowest Reynolds number feasible for testing, and they simulated approxi- 
mately the high-altitude portion of a Martian entry. The second major point 
of interest was the aerodynamic characteristics of the probe shape in the 
nose-forward attitude at high Mach number and Reynolds number simulating the 
low-altitude portion of the entry where the atmospheric measurements would be 
made. Drag and static and dynamic stability were determined from measurements 
of model attitude and position on shadowgraphs. The shadowgraphs were also 
used for flow visualization. The pitching -moment coefficients for the entire 
angle -of -attack range from 0° to l80° were determined by extrapolating experi- 
mental data for the nose-forward and nose-rearward attitudes and then improv- 
ing the extrapolations by an iterative technique based on matching a 
synthesized oscillatory motion to an observed free-flight motion. 
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mass of model, gm 
Mach number, dimensionless 
pitching velocity, radians/sec 
free-stream dynamic pressure, newtons/m 2 

Reynolds number based on maximum diameter and free-stream conditions, 
dimen s ionle s s 

flight time, sec 

velocity along flight path, m/sec 
initial velocity (at t = 0), m/sec 
distance along flight path, m 

axial distance from model base to center -of -gravity position, cm 

angle of attack (in the vertical plane), deg 

a max.j + a max.p 

average value of maximum-angle envelope, , deg 

(see sketch (a)) 

a min : ? + a minf 

average value of minimum-angle envelope, — g > deg 

(see sketch (a)) 

/ x 

a- 2 dx, deg 

s]a 2 + P 2 , deg 

angle of sideslip (in the horizontal plane), deg 
damping exponents in equation (l), mf 1 
wavelength of pitching oscillation, m/cycle 
free-stream air density, kg/m 3 

radius of gyration about a transverse axis through the center of 
gravity, cm 

dynamic-stability parameter, C-q - + (C^ + C^) 

dimensionless ^ " L 

roll parameter. r °l^- radians /m 
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root-mean-square angle of oscillation 
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w 1 2 rates of rotation of vectors which describe the model pitching motion 
3 in equation (l), radians/m 

( ) * first derivative with respect to x 

( )" second derivative with respect to x 

(’) first derivative with respect to time 


MODELS, TEST CONDITIONS, AND DATA REDUCTION 


Figure 1 is a detailed sketch of the model tested. The diameter d was 
0.508 cm for the models and 25*9 cm for the full-scale vehicle. To position 
the center of gravity of the models to match that of the full-scale vehicle, 
the models were constructed from two materials joined approximately at the 
center of gravity. The first models had aluminum afterbodies and sintered 
tungsten-iron -nickel composition forebodies. The two materials were held 
together by a steel machine screw. Initial flight tests revealed an undesir- 
able spallation of tungsten particles from the tip of the conical nose because 
of high heating rates; later models were constructed with tantalum- tungsten 
alloy forebodies to eliminate this problem. As seen in figure 1, the bases 
of the models had two cylindrical cavities which were also scaled from the 
full-sized vehicle. Figure 2 is a photograph of the models and the two types 
of sabots used - one for nose -forward tests and the other for nose-rearward 
tests. 

All models were gun-launched from deformable piston light-gas gun. Time- 
distance histories and attitude histories were recorded at spark shadowgraph 
stations along the flight path. Except for one flight in the pressurized 
ballistic range (which has 2k stations along a 62-meter length) at M ~ 10 
and R ~ 0.02x10 s , all flights in the nose -forward attitude were made in the 
prototype hypervelocity free-flight facility (which has 11 stations along a 
12.1-meter length) at Mach numbers of 10 and 26 and Reynolds numbers from 
0.55x10 s to 0.67x10 s . (One partial flight was made at M = 10 with 
R = 0.15x10 s .) The M = 10 tests were made in still air, and the M = 26 tests 
in a countercurrent airstream. 

Three nose-rearward flights were made in the pressurized ballistic range 
at Mach numbers from 7.0 to 9.2 and at Reynolds numbers from 0,009*10® to 
0.165x10 s . The physical characteristics of the models and the respective 
flight conditions for each test are summarized in table I. Tests were not 
conducted in gas compositions containing C0 2 , a more appropriate gas for a 
Martian atmosphere, because (l) the present investigation was not intended to 
be a comprehensive study of gas composition effects; (2) air as a medium is 
much simpler to use, especially at M = 26 where a countercurrent flow is 
necessary. 


k 



Stability Data -Reduction Methods 


Nose -forward tests .- Stability data in the nose-forward attitude were 
obtained from conventional analyses of the attitude histories of the models. 
Examples of the pitch-yaw motions experienced by the models are shown in 
figure 3 . Figure 3(a) shows a nearly planar motion and figure 3(b) shows a 
precessing elliptical type of motion. The numbered symbols show the angles of 
attack and sideslip measured from the shadowgraphs at each of the stations. 
(Motions shown represent flight segments of seven stations each.) The curves 
are fitted to the data points by a method which will be discussed below. A 
physical interpretation of this type of plot is to imagine oneself looking 
along the flight path of an approaching model and observing the path that 
would be traced by the nose of the model as it oscillated about its center of 
gravity and precessed about its flight path. The majority of the flights 
exhibited nearly planar motion with little precession. 

To obtain the static-stability parameter or quasilinear value of 
pitching -moment-curve slope, Cm™ > and - the dynamic -stability parameter 
(ref. 1) 1 

£ = C D “ + ( Cm q + C m ^)( ci / a ) 2 

from the attitude histories, the motion equation developed in reference 2 
(further discussed in refs. 3 and 4) was used. Specifically, the tricyclic 
equation 

o . • i (m+iwjjx 1 (rip - iw 2 )x , _ icp'x x 

3 + ia = kie' ' + k 2 e v ^ + k 3 e Y (1) 

was fitted by the method of least squares to the measured values of a and 3 
for each flight. The exponents 2 are related to the wavelength of 
oscillation by 

x = — 

■7“ i«*fe 

and. the exponents f)i 2 are related to | by 

| = + Tla 

pA/2m 

The curves in figure 3 represent best fits of equation (l) to the data. 

To identify the amplitude of oscillation of each flight for analysis and 
presentation of the data, both the minimum and the maximum angles of oscil- 
lation must be considered because the angle range through which the models 
oscillated differed for each flight. The amplitudes are defined in 
sketch (a). 
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Sketch (a) 

The method of data reduction using equation (l) is based on linear 
aerodynamics (thus a constant -wavelength. A); therefore, for bodies with 
highly nonlinear stability characteristics, these quasilinear measurements are 
further reduced by a method developed in reference 5, which takes into account 
the observed nonlinearities, and a curve of pitching moment versus angle of 
attack is obtained. A complete explanation of how the results of reference 5 
are applied to the present investigation is presented in the appendix. 

Where there was a sufficient amount of oscillatory motion in one flight 
and a marked difference in wavelength or drag was observed over the distance 
of the flight, the flight record was split into two separate overlapping 
flight segments for data reduction. Not all of the flights could be split in 
this fashion, but almost all of the M = 10 flights were. The exceptions are 
noted later. The typical motions seen in figure 3 ar e portions of two flights 
rather than the entire flights themselves. 

Nose-rearward tests . - Since the model was presumed initially to be 
statically unstable in the nose -rearward attitude, it was decided that the 
pitching moment would have to be deduced directly from the pitch acceleration 
of the model as it rotated away from the unstable attitude and toward the sta- 
ble position of a = 0°. Because of restrictions on gun placement with 
respect to the shadowgraph stations and unavoidable sabot interference upon 
launch, the initial pitching motion of the model from the a = 180° position 
in the gun could not be measured, so the test was designed to obtain as many 
shadowgraphs as possible the second time the model was in the vicinity of 
a = 180° (one-half cycle of motion from the initial launch position). A well- 
defined angle of attack versus flight -distance curve was obtained with a least 
squares fit of a polynomial in x to the data in the vicinity of a = 180° 
(either at the first observable amplitude peak of an oscillatory motion about 
a = 0° or at a = l80° if the model passed through a = l80° in a slow tum- 
bling motion produced by unsymmetric sabot separation). The equation of 
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motion governing the behavior of the model, if aerodynamic damping is 
neglected and velocity loss is small is 


where 


la" 


2 


'm 


a' 


dfa 

dx 5 


Wear peak amplitude, a* is small, and the aerodynamic damping and drag 
effects (which are proportional to a*) can quite properly be neglected. If 
the variation, a = f(x) (where x is flight distance) can be determined to 
represent the motion, it can be differentiated twice to obtain 


a" = f"(x) = C 


m 


or 


Cm 


- Si 


We therefore have equations for a and C m as polynomials in x. These equa- 
tions then can be used to produce a curve of C m versus a for the range of a, 
being considered. 


RESULTS AND DISCUSSION 


Plow Fields 

Pertinent features of the' flow fields around the models will be described 
first for later reference in a discussion of the stability data. 

Shadowgraphs of a model in the nose -forward attitude at M = 10 and 
R = 0.6xl0 6 at angles of attack from 0° to 2° are presented in figure 4(a). 

The shadowgraphs clearly show flow separation over the cylindrical part of the 
body and part of the flare. At a = 0° the point of separation appears to be 
slightly behind the cone-cylinder junction. As the angle of attack increases 
to 2°, the flow over the cylinder and flare gradually becomes essentially 
attached on the windward side and fully separated on the leeward side. 

Another important characteristic is that the flow over the entire body is lam- 
inar, as can be verified by examination of the wake detail in the original 
pictures. 

The observed bulge in the bow-shock wave is caused by slight blunting of 
the conical tip. (The model tip was blunted slightly initially to conform to 
full-scale geometry, and became slightly blunter in flight because of erosion 
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of the heated nose. Erosion is greater in the tests at M = 2 6 . These 
shadowgraphs are all from the same flight of a model with a tantalum-tungsten 
nose. There is very little change in tip geometry over the length of the. 
flight. For comparison, figure 4(h) shows shadowgraphs of a model at M = 10, 
R = 0.6xl0 6 with a sintered tungsten-iron-nickel composition nose. The first 
shadowgraph was taken early in the flight and shows laminar separated flow 
over the cylinder, as does figure 4(a). The second shadowgraph was taken at 
the last station in the flight and shows that the nose material has begun to 
melt, creating what appears to be turbulent flow, or possibly laminar flow 
containing liquid iron and solid tungsten particles. The nose has been sig- 
nificantly blunted.. The tests involving these models will be pointed out 
later in discussing the data and will be referred to as those with "turbulent" 
flow. 

Figure 5 shows a series of four shadowgraphs of a model at M = 26 and 
R = 0.6x10 s ranging in angle of attack from 8° to l4°. The major point illus- 
trated is the movement of the bow-shock wave toward the windward side of the 
flare as the angle of attack is increased. At a near 12°, the bow shock 
appears to impinge on the flare. The streak in the shadowgraph for a = l4.2° 
is a streak in the test section window and not a vortex streak as might be 
suggested. 

Figure 6 shows a group of shadowgraphs of a model at M =26, R = 0.6xl0 6 , 
and low angles of attack. Although the report reproductions may not be suf- 
ficiently clear, the original negatives show evidence of separated flow, 
similar to that observed at M =10. This similarity in observed flow fields 
suggests that over this Mach number range, at low angles of attack, the flow 
characteristics are insensitive to Mach number. This provides some 
justification for analyzing the M = 10 and 26 data together. 


Drag 

Drag coefficients were computed for each flight by least-squares fitting 
the following equation to the time, t, and flight distance, x, measurements 
recorded at each data station 


* - - « <« f - 6 > 

where k = (pA/2m). The drag coefficient, Cp, and initial velocity, Vq, are 
the coefficients determined from the fit. Figure 7 shows the drag coefficient 
plotted against a rms , the root-mean-square angle of attack, for all flights. 
Each of the flights at M = 10, except the two with models that had spalling 
nose material and the one conducted at low Reynolds number, was split into two 
flight segments for data-reduction purposes because of the large variation of 
Cp with a rms . The M = 26 flights did not experience enough cycles of motion 
to be split. The data point at low Reynolds number was obtained from a flight 
at low air density with very limited time-distance information; the degree of 
uncertainty in this point is illustrated by the scatter bar. A curve has been 
faired through the data points in figure 7* For comparison, a least-squares 
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fit of a quadratic equation to all the data except the three filled points is 
shown also. A quadratic fit, commonly used to represent drag data for linear 
systems, is seen to he a rather poor representation of the Cp data at the 
larger amplitudes . Also shown is Newtonian theory which overpredicts the drag 
at a = 0° by about 60 percent. 

To obtain a possibly more precise extrapolated value for at a = 0° 

( C j) Q ) , the data points below cx rms = I 3 0 were fitted in three different ways. 
The results are shown in figure 8 . Again the filled symbols are not included 
in the least-squares fits. The Mach number 10 data, fitted alone, appear to 
have a slightly lower Cp Q and a faster rise with a rms than the Mach number 
26 data, fitted alone, or the all -data curve; however, caution should be used 
in interpreting this difference as a Mach number effect. A more striking fea- 
ture is that the two data points below a = 2 ° are low with respect to the 
quadratic fits through the data. This is not surprising considering that the 
flow separates over part of the flare and cylinder at a = 0 ° and becomes 
fully attached on the windward meridian at a = 2°. Shown also are two calcu- 
lated values for Cp 0 - one obtained from embedded Newtonian theory (ref. 7) 
and one based on a separated-flow calculation. The separated-flow estimate is 
based on conical flow theory for the cone, Prandtl -Meyer expansion at the 
cone -cylinder junction, and oblique shock relationships through the shock 
waves where flow separation begins and where the flow reattaches on the flare. 
The separated-flow region is considered to be a solid surface or dead-air 
region with identical pressure to that behind the shock wave at the separation 
point on the cylinder (see sketch (b)). The geometric locations of separation 



and reattachment points were measured from a shadowgraph of the model at 
a = 0° and M = 10 (see fig. 4(a)). Skin-friction drag on the cone was 
included in the separated-flow drag estimate. 


Static Stability - Nose -Forward Attitude 

The quasilinear data for static stability, Cn^, both for M = 10 and 26 
were considered as a group in deducing the pitching -moment-coefficient curve 
since there was no indication of Mach number effects either in the computed 
results or in the flow fields observed in the shadowgraphs. The quasilinear 
results are recorded in table I. All flights at M = 10 were divided as 
mentioned earlier, except four that either did not have enough cycles, or had 
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such small amplitudes that no improvement would result. None of the flights 
at M = 26 contained enough cycles to permit division into parts. The method 
of reference 5 was used to obtain a pitching -moment curve (explained in detail 
in appendix) and a plot of C m versus a is shown in figure 9* Because of the 
high degree of nonlinearity over a wide angle-of -attack range , two separate 
polynomials were required to represent the pitching -moment characteristics 
adequately and the region of overlap is shown as the small area between the 
curves. 

The static stability (represented by the local slope of the pitching- 
moment curve) is highest in the region of a = 0° and then drops off rapidly 
to a minimum at a - 6° , where the stability begins to increase again and 
reaches a secondary maximum value in the vicinity of a = 12° to 14°. This 
region of high stability at a near 12° correlates quite well with the previ- 
ously discussed bow-shock impingement on the flare. The high stability at 
a = 0° can be attributed to the separated laminar boundary layer over the 
cylinder and flare observed at low angles of attack. Reference 8 shows that 
laminar separation over flared bodies of revolution results in increased 
normal-force -curve slope and large rearward movements of centers of pressure. 
Similar results are presented in reference 9. 

Newtonian theory is also shown for comparison. It considerably under- 
estimates both the initial slope and the magnitude of the pitching -moment 
curve . 


Static Stability - Nose-Rearward Attitude 

Two flights were made for the purpose of determining pitching moments at 
angles of attack near 180°. Both flights were made at low-air densities to 
achieve as long a wavelength of oscillation as possible and thereby obtain a 
maximum number of shadowgraphs of the models in the attitude of interest. The 
low densities also resulted in Reynolds numbers that were in the range of 
interest for full-scale nose-rearward flight. It should be noted that both of 
the nose -rearward flights were nearly planar in nature with no measurable roll. 
The first flight was at M = 7.82 and R = 0.0108xl0 6 . The resulting motion is 
shown in figure 10 (a), where the variation of angle of attack with flight dis- 
tance is plotted. As seen in the figure, the model oscillated to a^x = l65° 
after a half cycle of motion from its attitude of l80° in the gun muzzle. The 
filled symbols were used to determine the equation of a = f(x). The actual 
maximum value of a and its location along the flight path were found first 
by least squares fitting the ten filled symbols with a fourth -degree polynomial 
in x and then setting the first derivative of the polynomial equal to zero. 
Two segments of the curve, one on each side of the maximum point, were then 
fitted by least squares with fourth-degree polynomials in x and differenti- 
ated twice to get a” = f Tt (x). The pitching moments derived from these two 
segments were then plotted versus angle of attack and considered to be accu- 
rate over the range of angles of attack where the respective C m values from 
the segments agreed reasonably well. They agreed within 2 percent from 
a = 165 0 to 130 o , so this was considered to be the region of confidence. 
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The second flight was made at the lowest density possible for the 
facility at M = 9.2 and a resulting R = O.OO 89 XIO 6 . The motion experienced 
during this flight is seen in figure 10(b) . The model rotated one complete 
revolution and passed through a = l80° with a small a T . The filled symbols 
were used to obtain the desired pitching moment. The first objective was to 
determine the value of a 1 at a = l 80 ° and the corresponding value of x. 
This was done by least squares fitting a third-degree polynomial in x to the 
three points directly below and the two points directly above a = l80°. Once 
a 1 was determined, the procedure was essentially the same as for the previous 
test, except that the boundary value of a* was not equal to zero. The 
derived pitching moments from the two curves agreed very well from a, = l 80 ° 
to 165 °. The pitching -moment curve shown in figure 11, extending from 
a = l80° to I 3 O 0 , was derived from those two flights. The Newtonian theory 
prediction, also shown, compares favorably with the experimental result in 
this angle -of -attack range. 


Construction of Entire Pitching -Moment Curve 

With C m determined experimentally over ranges of angles of attack near 
0° and l80°, one could construct a reasonable C^ curve for the complete 
angle -of -attack range by choosing a theoretical value in the vicinity of 
a = 90° to 100° (Newtonian value, for instance) and by extrapolating the 
experimental curves through the chosen point. However, it is possible to 
obtain a more precise definition of the pitching -moment curve near a = 100° 
and, consequently, of the entire C m curve, by an iterative technique in the 
following manner. An additional nose -rearward launch was made at M = 7 and 
R = 0.165X10 6 in which the model oscillated through six cycles of motion, and 
damped from a maximum amplitude of l45° to 85°. (ideally one would like to 
launch a model at an initial angle of a = 9 O 0 to 100° and analyze the motion 
just as was done for the a = l 80 ° region, but model length prevented placing 
the model in the gun in this position. ) The Reynolds number was higher than 
in the other two flights when the models were launched nose rearward because 
of the higher free -stream density necessary to produce more cycles of motion 
and because of the resultant damping to the amplitude range of interest. A 
synthetic motion, generated from numerical integration of the equations of 
free -flight motion, was matched to the observed motion using as inputs the 
flight conditions of the observed test (i.e., free-stream density and velocity 
history), the measured model characteristics, and the pitching -moment curves 
derived experimentally for the nose-forward and nose-rearward attitudes. (See 
ref. 10 for a complete derivation of the equations of motion and a description 
of the motion-synthesis program. ) Starting values of pitching moment for 
angles of attack between 40° and I 3 O 0 are determined by smooth extrapolation 
of the experimental data for nose-forward and nose-rearward attitudes to a 
value predicted by Newtonian theory for a* = 100° (see fig. 12). This esti- 
mated value and the resulting extrapolations were then altered by trial and 
error until the synthesized motion matched the observed motion. The con- 
structed C m curve was altered simply by raising or lowering the maximum 
point near a = 100 ° and maintaining a smooth interpolation from each of the 
previously established curves. Figure 12 shows the final constructed curve 
for C m versus a, as well as the Newtonian curve. Figure 13 shows the 
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synthetic motion and its relation to the data points from the observed motion. 
The value shown for CMQ, (g - Cp), a modified dynamic -stability parameter 
(assumed constant for the entire angle range), is the final iterated value 
which in conjunction with the selected pitching -moment curve produces the best 
match to the observed amplitude change. 

Figure l4 compares the observed data points and a motion generated with 
Newtonian pitching -moment coefficient. The match is seen to be poor, as the 
resulting frequency is somewhat low. 


Dynamic Stability, Nose Forward 

The dynamic-stability parameter £ is plotted in figure 15 versus a m , 
the average maximum amplitude in pitch, for each flight or flight segment. 

The model shows good dynamic stability throughout the angle range tested. 
Although several of the free -flight motions (including all flights at M = 2 6 ) 
are not adequate for defining dynamic stability accurately, the trend is for 
dynamic stability to decrease with increasing angle of attack. The models 
with "turbulent" flow (filled symbols) are slightly less stable dynamically 
than the models with laminar flow (open symbols). The large values of |, 
obtained for low angles of attack, may be strongly influenced by the nonlinear 
static stability, which increases rapidly as the angle of attack decreases, 
or they may simply indicate poor definition of amplitude in the small angle 
range where errors of tenths of a degree become significant. 


CONCLUSIONS 


Drag and stability characteristics and flow-field observations were 
obtained from free-flight tests of models of the slender flare-stabilized cone 
cylinder proposed by Goddard Space Flight Center for a Martian atmosphere 
probe. The results presented herein, for tests conducted in both a nose- 
forward and nose -rearward initial attitude, lead to the following conclusions: 

1. Static stability decreases rapidly with increasing angle of attack 
until the bow-shock wave approaches the windward side of the flare. Further 
increase in angle of attack produces a notable increase in stability. The 
reduction in stability with increasing angle of attack at low angles was due, 
at least in part, to the observed laminar -boundary-layer separation over the 
cylinder and flare at low angles. 
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2. There are no significant differences in drag, stability, or flow- 
field characteristics between Mach numbers of 10 and 2 6 at a constant Reynolds 
number of 600,000; no effect of Reynolds number on drag or stability was 
measured at Reynolds numbers from 23,000 to 600,000 at a Mach number of 10. 

3. The configuration is statically unstable in the nose -rearward 
attitude. The only stable trim point is at a = 0° (nose forward). 

4. The configuration is dynamically stable at all angles of attack. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Jan. 20, 1967 
124 -07 -02 -13 -00 -21 
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APPENDIX A 


METHOD OF OBTAINING NONLINEAR PITCHING -MOMENT 
CURVES FROM QUASILINEAR DATA 


For a pitching -moment equation of the form of a polynomial. 


■C = P o a + P^jx 2 + P 2 a 3 + P3a 4 + P 4 a 5 + P5CL 6 + P 6 a 7 


(Al) 


the equation derived in reference 5 for the quasilinear value of the pitching- 
moment-curve slope can be written as 




= P D a + | Pib + 1 P 2 c + | P 3 d + I P 4 e + | P 5 f + + 
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'm min 


d = 


a 2 . a 2 \ 7/2 
a m + a min 


(a 2 - a 2 . Y 
v m min' 


e = 


7a 4 + 10a 2 a 2 . + 7a 4 . 

; m m min 1 min 


a® + a. 9 


f = 


m min _ / m 

2 V 2 


a 2 + a 2. \9/2 


(a 2 - a 2 . ) 

' m mm -n * 


min' 


g = 3a 6 + 5a 4 a 2 . + 5a 2 cc 4 . + 3a 6 . 

0 m m min m min J min 


A computer program was developed for applying these equations. It uses 
the method of least squares, considers all possible polynomials with from two 
to five of the seven terms shown in equation (Al), and selects the one poly- 
nomial from each of these four classes of polynomials (i.e., two, three, four. 
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or five terms) that produces the best fit to the , expand data. 

The class of polynomials judged to be the most representative of the data is 
then chosen. (The best five -term polynomials -will always be a better fit than 
the best four -term polynomials, for example, but may not be a reasonable 
representation of the data. For example, a five term polynomial could be made 
to fit five data points exactly, but the resulting curve might have such large 
variations between points that it would be completely useless. A three- or 
four-term polynomial, on the other hand, might prevent these large variations 
and be a more reasonable curve even though it can not pass through each point 
exactly.) With the coefficients P Q , P-l, P 2 , etc., determined, the pitching- 
moment coefficient is expressed by equation (Al), where some values of P 
will be omitted, depending on the number and degree of terms. 

Figure 1 6 illustrates the application of this procedure to the present 
data. From the variation of the apparent static stability with angle of 
attack, it is obvious that the static stability is not linear. It can be seen 
from equation (A2) that if n / 0 and the pitching -moment slope is not 
linear, there is no exact form for plotting directly on a two- 

dimensional plot (three dimensions are required). For most test motions, how- 
ever, is much smaller than a m and an approximate representation of the 

data (such as fig. 1 6 ) can be shown by plotting ve rsils a m equa- 

tion (Al), where cxjnin is set equal to zero. Because the variation of the 
data with a m was not simple, two separate polynomials were necessary to 
represent the data adequately. These two curves in their approximate form 
(approximate because ami n bas been set equal to zero) overlap in the vicin- 
ity of 0 ^ = 15°. These two polynomials produce the two polynomials 
describing Cm = C m(°0 shown in figure 9« 
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TABLE I.- MODEL CHARACTERISTICS AND TEST CONDITIONS 


Test 

(a) 

M 

RXlO -6 

m/sec 

kg/m 3 

a m^ 

deg 

a min> 

deg 

a rms > 
deg 

c D 

1 

0 

£ 

“1 

m, 

gm 

<L 

cm 

_C£ 

d 

IxlO 8 , 

kg-nf 2 

(i)‘ 

Facility 

0>) 

c l463 

9.87 

0.606 

3378.7 

O.636O 



4.30 

0.144 



1.1407 

0.5098 

2.369 

3-38 

0.876 

PHFFF 

0 1463-1 

9.89 

.607 

3385.4 

.6360 

6.49 

0.23 



1.847 

7-3 

1.1407 

.5098 

2.369 

3-38 

.876 

PHFFF 

c l463-2 

9.85 

.605 

3372.0 

.6360 

5.98 

.24 



1.950 

3.6 

1.1407 

.5098 

2.369 

3.38 

.876 

PHFFF 

°1594 

10.38 

.632 

3554.0 

.6329 

2.12 

•99 

1.63 

.102 

2.627 

20-9 

1.1453 

.5080 

2.400 

3-46 

.855 

PHFFF 

885 

10.42 

•0231 

3513-1 

.02288 

31.34 

1-93 

22.96 

.850 

1.793 

20.8 

1.1267 

•5075 

2.376 

3.28 

.884 

PBR 

1521 

9-76 

.148 

3362.5 

.1578 

16.08 

.98 

10.66 


1-597 


1.2052 

• 5085 

2.424 

3*44 

.904 

PHFFF 

1542-1 

9.97 

.601 

3432.3 

.6272 

7.00 

.20 

4.98 

.136 

1.415 

28.2 

1.1647 

•5093 

2.389 

3*44 

.876 

PHFFF 

1542-2 

9-94 

.599 

3421.7 

.6272 

4.61 

.28 

3.04 

.111 

1.866 

55-2 

1.1647 

•5093 

2.389 

3.44 

.876 

PHFFF 

1668 

10.22 

.603 

3544.2 

.6194 

2.35 

.35 



3.090 

73-2 

1.1607 

.5070 

2.387 

3.34 

.894 

PHFFF 

1668-1 

10.23 

.604 

3547.9 

.6194 



1.79 

.093 



1.1607 

.5070 

2.387 

3-34 

.894 

PHFFF 

1668-2 

10.21 

.603 

3540.2 

.6194 



.88 

.079 



1.1607 

.5070 

2.387 

3-34 

.894 

PHFFF 

1669-1 

10.45 

.632 

3597-8 

.6282 

5.80 

1.05 

4.23 

.122 

1.779 


I.1808 

.5105 

2.378 

3-37 

.914 

PHFFF 

1669-2 

10.42 

.631 

3588.4 

.6282 

5.06 

.74 

3.60 

•113 

1.936 

12.2 

1.1808 

•5105 

2.378 

3.37 

.914 

PHFFF 

1671-1 

10.31 

• 573 

3584.1 

.5832 

5.32 

2.32 

4.08 

.130 

1.837 

3.8 

1.2045 

.5065 

2.428 

3.43 

.901 

PHFFF 

1671-2 

10.29 

• 571 

3573-2 

.5832 

4.77 

1-93 

3-63 

.128 

1.852 

6.4 

1.2045 

• 5065 

2.428 

3.43 

.901 

PHFFF 

1672-1 

10.39 

.563 

3657.0 

.5690 

9.22 

.21 

6.21 

.158 

1.570 

8.3 

1.1698 

•5098 

2.387 

3-40 

.896 

PHFFF 

1672-2 

10-37 

.561 

3648.1 

.5690 

8.79 

.48 

6.36 

.156 

1.564 

8.4 

1.1698 

.5098 

2.387 

3-40 

.896 

PHFFF 

1673-1 

IO.38 

• 578 

3516.2 

•5934 

46.71 

5.10 

32.69 

1.235 

1-555 

7-7 

1.1757 

.5108 

2.368 

3.43 

.894 

PHFFF 

1673-2 

9.87 

.564 

3421.7 

•5934 

38.75 

3.80 

26.72 

.990 

1.688 

9.8 

1.1757 

.5108 

2.368 

3.43 

.894 

PHFFF 

1547 

26.7 

• 552 

5197-0 

.1382 

9.82 

2.15 

6.87 

.154 

1.007 

17.5 

1.1906 

.5080 

2.4l8 

3-39 

.906 

PHFFF 

1548 

26.7 

.614 

5266.0 

.1534 

12.14 

2.50 

8.58 

.186 

1.430 

10.1 

1.1888 

.5126 

2-379 

3.46 

.903 

PHFFF 

1549 

25.6 

.613 

5100.0 

.1621 

15.35 

.11 

10.46 

.235 

1-563 

16.2 

1.1718 

.5105 

2.385 

3-35 

.911 

PHFFF 

1550 

26.2 

.626 

5225.O 

.1618 

15.69 

1.93 

11.21 

.243 

1-398 

24.9 

1.1650 

•5093 

2.371 

3-31 

.913 

PHFFF 

1551 

26.0 

.625 

5186.0 

.1629 

18.84 

1.04 

12.90 

.305 

1-750 

7.2 

1.1795 

• 5105 

2.380 

3.43 

.896 

PHFFF 

1552 

26.2 

.666 

5225.0 

.1720 

13.80 

1.04 

9.47 

.220 

1.360 

30.4 

1.1802 

.5118 

2.387 

3.43 

.900 

PHFFF 

d 842 

7.03 

.165 

2376.5 

.2422 







1.1436 

.5088 

2-377 

3.34 

.887 

PBR 

d 974 

7.83 

.0107 

2652.0 

.01422 







1.1321 

• 5095 

2.371 

3.36 

.874 

PBR 

d 975 

9.22 

.00889 

3144.9 

.009947 







1.1570 

• 5105 

2.383 

3-39 

.890 

PBR 


^Ihe numerals 1 and 2 following test numbers indicate first or second segment of a complete flight that has been split into 
two separate parts for analysis. (Explanation in section on data-reduction method.) Test 1463 was split for stability data but 
not for drag. Test 1668 was split for drag but not for stability. 

^PHFFF - Prototype Hypervelocity Free -Flight Facility; PBR - Pressurized Ballistic Range. 
c Models with "turbulent" flow. 

^-Models launched in nose -rearward position. 

Note: Blanks indicate parameter is not applicable for that flight or flight segment. 
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Figure 2.- Models and sabots. 


A-37592 







a = 2.0° (p = -0.6°), M = 9.57, R = 0.600X10 S A- 36957-4 

(a) Test 1542* tantalum -tungsten alloy nose. 

Figure 4.- Shadowgraphs of model in flight at M = 10. 
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a = 0.8° (p = -5-0°), M = 9-92, R = 0.609X10 6 


A-34745 



a = -0.6° (p = 0 .0°) , M = 9.82, R = 0.603X10® 

(b) Test 1463, sintered tungsten-iron-nickel nose. 
Figure 4. - Concluded. 
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' ■■.■*:* . > 
t ' '*» " : 


-- •*: 7%^-j ;/ i: 

^3pyA;iL;;ljijSlI:' : i 


a 


= -10.6 (|3 = -1.0 ), M = 26.2, R = O.67OXIO 6 


A-36957-6 



Figure 5.- Shadowgraphs of model at M = 26 at various angles of attack. 
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Figure 8.- Variation of drag coefficient with a^g and estimate of Cj) 0 . 


0.01771 a + 0.006073 a 2 - 0.0002626 a 3 + 0.000004766 a 4 - 0.0000000323 a 



-C m =0.06923 a -0.01 144 a 2 +0.0008808 a 3 -0.00001990 a 4 











300 



t, sec 

Figure 13.- Comparison of observed data points (test 842) and motion synthesized with experimental 

C m curve. 
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Figure 14. - Comparison of observed data points (test 842) and motion synthesized with Newtonian 

C m curve. 
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Figure 15.- Variation of dynamic -stability parameter with average maximum 

amplitude of oscillation. 
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Figure l6.- Variation of quasilinear pitching -moment curve slope with average maximum amplitude of 

oscillation. 
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